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In this paper we construct small amplitude periodic internal waves traveling at the boundary region 
' between two rotational and homogeneous fluids with different densities. Within a period, the waves 

we obtain have the property that the gradient of the stream function associated to the fluid beneath 
' the interface vanishes, on the wave surface, at exactly two points. Furthermore, there exists a critical 

\ layer which is bounded from above by the wave profile. Besides, we prove, without excluding the 

^ . presence of stagnation points, that if the vorticity function associated to each fluid in part is real- 

^ ■ analytic, bounded, and non-increasing, then capillary-gravity steady internal waves are a priori 
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on the vorticity function. 
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1. Introduction 



In this paper we consider two-dimensional internal periodic waves traveling at the interface 
between two layers of immiscible fluids with different densities, under the rigid lid assump- 
tion. In our context, the fluids have constant vorticity and we construct internal travel- 
J> \ ing waves with a critical layer and stagnation points in both gravity and capillary-gravity 

■ regimes. 

■ A critical layer is a region of fluid consisting entirely of closed streamlines, and stagnation 
. points are fluid particles traveling horizontally with the same speed as the wave. Flows with 

stagnation points are known [5] to be relevant for the description of background states for 
tsunamis. Concerning waves traveling over a homogeneous fluid, it was observed in [SU^ that 
the strong elliptic maximum principle rules out the existence of smooth irrotational waves 
with stagnation points or critical layers. On the other hand, it was shown in |31H33) that 
k> . there exist extreme Stokes waves which are Lipschitz continuous and have stagnation points 

}^ ' and a sharp corner at the crest. For the rotational case the picture is different. Existence of 

^ ■ exact periodic traveling gravity waves with general vorticity has been established first in [6] 

by means of bifurcation and degree theory, and in [8] in the case of waves with bounded 
and discontinuous vorticity. By construction, the waves in [Bl[8] do not possess stagnation 
points and critical layers. Linear gravity water waves with stagnation points which travel on 
currents with constant vorticity were studied in |15| , and the exact picture of the waves has 
been obtained first in |34j . and later on in |9] by employing complex methods. These waves 
possess at most a critical layer which is located within the fluid body. 

One of the factors which can determine the presence of critical layers is the vertical 
stratiflcation of the fluid. It was recently shown in |19| that continuously stratified gravity 
waves with constant vorticity may possess two critical layers and the qualitative picture of 
the streamlines may be different than that for homogeneous flows [2l[Sll- Another factor is 
the presence of an affine vorticity distribution, situation analysed in [13pi4). when the waves 
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may possess arbitrarily many critical layers. For a survey on traveling waves with critical 
layers we refer to the article |16) . 

Internal waves are usually created by the presence of two different layers of water com- 
bined with a certain configuration of relief and current. They form where the water above 
and below the interface is either moving in opposite directions or in the same direction at 
different speeds. A well-known example are the internal waves which move from the Atlantic 
Ocean to the Mediterranean Sea, at the east of the Strait of Gibraltar. In this case the two 
layers correspond to different salinities, whereas the current is caused by the tide passing 
through the strait. The mathematical theory of waves at the interface between two layers 
of immiscible fluid of different densities has attracted a lot of interest and we refer to the 
survey article of Helfrich and Melville |20| which provides a good overview on steady internal 
solitary waves in such systems. 

In this paper we construct periodic steady internal waves between two layers of homoge- 
neous fluids with different densities in the gravity and capillary- gravity regimes. Particularly, 
we find solutions which possess exactly one critical layer which is bounded from above by the 
wave profile. Moreover, the gradient of the stream function associated to the fluid beneath 
the interface vanishes at exactly two points on the interface, meaning that, close to these 
points, fluid particles located below the interface move almost horizontally with velocity 
approaching the wave speed. The fluid above the interface does not sense the presence 
of these stagnation points. In contrast to the situation in |19) . where waves with critical 
layers were constructed in an unstable regime, the solutions we find are stably stratified. 
While in [31|I33| the wave surface is only Lipschitz continuous in a neighbourhood of the 
stagnation point, our solutions are real-analytic. In the pure gravity case this follows by 
using a regularity result for free boundary value problems, cf. [26], whereas when capillary 
plays a role we provide a new idea based on parabolic theory. It is worthwhile to mention 
that our method may be applied to prove the a priori analyticity of water waves traveling 
over an homogeneous fluid, when capillary effects are incorporated, provided the vorticity 
function is bounded, analytic, and non-increasing even if stagnation points exist (situation 
complementary to that analysed in |2ip23)). 

We determine, by using higher order expansions and elliptic maximum principles, a 
precise picture of the streamlines in the fluids which could be used to describe the particle 
trajectories. Due to the analyticity of the interface, one could choose also linear theory, 
similarly as in |10 1 [T5 l [25 l l29). to obtain an approximative picture of the exact particle paths. 

The outline of the paper is as follows: in Section [2] we present the mathematical model, 
re-express the problem in terms of the stream functions and provide the main regularity 
result Theorem 12.11 In Section [3] we write the problem as a nonlocal equation for the wave 
profile, and use bifurcation theory to prove the existence statement Theorem 13.11 as well as 
higher order expansion for the bifurcation curves. In the last section, based on Theorem 
12. 1| we illustrate the precise picture of the streamlines for some of the waves obtained in 
Theorem 12.11 



2. The governing equations and a regularity result 

In this section we present the governing equations for two-dimensional internal water waves 
travelling at the boundary region between two rotational fluids with different densities under 
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the rigid lid condition at the top. The bottom of the ocean is assumed to be flat and we 
denote with S the unit circle, i.e. S stands for R/27rZ. The fluid domain Q := S x ( — 1, 1) 
contains two fluids layers separated by a sharp interface y = r]{t,x), the wave profile, which 
defines the two subsets 

0^ := {{x,y) : X £ E and — 1 < y < ri{t,x)}, il* := {{x,y) : x G S and r]{t,x) < y < 1}. 

The domain il^ contains a Newtonian fiuid with constant density p, velocity field [u, v), and 
pressure P, and we denote by J), {u,v), and P the density, velocity, and pressure of the fluid 
located at the top. The line y = —1 is the impermeable bottom of the ocean and y = 1 is 
the rigid top of the two-fluid system. Assuming that both fluids are inviscid, the dynamics 
can be described by Euler's equations (see |27) for a justification of the inviscid fiow). Being 
interested in traveling internal waves, we presuppose that there exists a positive constant c, 
the wave speed, such that ri{t,x) = rj[x — ct), 

{u,v,P){t,x,y) = {u,v,P){x - ct,y), {u,v, P){t,x,y) = {u,v,P){x - ct,y), 

and formulate the problem in a frame moving with the wave. The equations of motion are 
the steady state two-dimensional Euler equations 

' {u- c)ux + vuy = -Px/p in il^, ( {u - c)ux +vuy = -Px/p in ^J* , 
< {u - c)vx + vvy = -Py/p - g in n'^i, < {u - c)v:^ + vvy = -Py/p - g in ^l^^, 
Ux + Vy = in n^, [ u^-\-Vy = in Q^, 

(2.1a) 

subjected to the boundary conditions, see e.g. [2l[6], 

. P = P- an" /(I + 77'2)3/2 ony = ri(x), 
V = U on y = —1, ^ I , ' ; : n ^ 

and < V = [u — c)rj on y = ri[x), (2.1bj 



V = on y = I, 



V = {u — c)r]' on y = ri{x), 



with 0" > being the surface tension coefficient and g the constant of gravity. To be more 
precise, we are interested in finding solutions of problem ()2.ip within the class 



7?gC2+"(S), {n,v,P)e (ci+°(05^))', (n,U,P)e (Ci+"(Q*))' (2.2) 

for some a € (0, 1). Similarly as in [6], we reformulate the problem ()2.ip by introducing the 
stream functions : Q,'^ ^ M. and ip : defined by the relations 

{u{x , s) — c) ds and ■il;{x,y) := / {u{x , s) — c) ds . 

1 Jri(x) 

Here m is a positive constant fixed such that ip = on y = rj[x). Indeed, taking into account 
that V'i/' = {—V, u — c) and Vip = {—v,u — c) it follows, by using the chain rule and ()2.1bp . 
that tp is constant on y = r}{x) and we can make this choice. Similarly, ip is constant on 
the rigid lid y = 1, and we let m := ip{Q, 1) be this constant. The above properties show 
that the streamlines coincide with the level curves of the stream functions. Both fiuids being 
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rotational, we introduce their vorticity by 

LO := Uy — Vx = Alp in f]^ and oJ ■.= Uy — Vx = Ai/j in $7^. 
Next, we assume that 

u - c < in ^ and u - c < in ^, (2.3) 

and exclude so the presence of stagnation points in the fluids. However, this assumption 
()2.3p will be dropped later on, after expressing the problem ()2.ip in terms of the stream 
functions. This will allow us to find solutions of system (|2.ip which do possess stagnation 
points. 

Condition (j2.3p guarantees [6l[30] that the vorticity is a single- valued function of the 
stream function, that is, there exist 7 G C"([— m, 0]) and 7 G C^([0,— m]), the vorticity 
functions, with 

uj{x,y) = j{-i;{x,y)), (x,y)Gri* and u]{x,y) = j{-'il){x,y)), {x,y) € n''^. 
Since, by Bernoulli's principle, the quantities 

(ti-c)2+t;2 P , J {u-cf+v^ P f^_, 

\ Vgy - / i{-s) ds and V — + gy - / i{-s) ds 

^ P Jo ^ P Jo 

are constant in fi* and f]^, respectively, we find, by restricting to y = r](x) that 

1" 52/ H = const. and \- QV + — = const. 

2 p 2 p 

To conclude, we observe that ip and ip are solutions of the semilinear Dirichlet problems 

C AV^ = 7(-V') in fi^, r A^ = 7(-V^)inn*, 

ip = ^ on y = r/(x), and < -0 = on y = r/(x), (2.4a) 
ip = m on y = —1, tp = fn on y = 1, 

respectively, and they are coupled by the following equation 

p\Vlp\'^ - p\V^\^ + 2g{p- p)y + 2ari"/ {1+7]'^ f/^ = Q on y = r]{x) (2.4b) 

for some constant Q G M. If we integrate this last equation ()2.4bp over the unit circle, we 
obtain, by taking into account that 'r]"/{l + _ (^r]'/{l + r/'^)^/^)', that 

Q = [ ip\Vi^Wx, r]{x)) - p|VVp(x, r]{x)) + 2gip - p)r?(x)) dx. (2.4c) 
Js 

The integral over S is normalised such that J^ldx = 1. It is not difficult to check that 
given two vorticity functions 7 and 7 and a classical solution {rj^ip^ip) of (j2.4p . then we can 
associate to this solution a unique solution of problem (|2.ip . even if the condition ()2.3p is 
not satisfied. This formulation of problem (|2.ip is much more convenient because we deal 
with two coupled semilinear elliptic problems. If 7 and 7 are such that we may solve (|2.4ap 
for given 77, which is definitely the case when the fiuids have constant vorticity, then (|2.4p 
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reduces to the problem of determining rj. As in |34) . we remark that the solutions of (j2.4p 
are solutions of (j2.ip for any value c G (0, oo) of the wave speed. 

As a first result, we use in the pure gravity case o" = a theorem on the regularity of free 
interfaces from |26| . which was employed in [Hl21 |[22tl30] to establish analyticity of traveling 
water waves without stagnation points on the profile in different wave regimes, to prove that 
internal traveling waves with real- analytic vorticity functions are a priori real- analytic. When 
we incorporate surface tension effects, we provide a new idea, which is based on parabolic 
theory, and show the real-analyticity without excluding the existence of stagnation points. 

Theorem 2.1. Let a € [0, oo) and assume that 7 and are constant functions. Given a 
solution of (j2.4p within the class (|2.2p . we assume, when a = 0, that 

IVV'P + |VV^P > ony = 7]{x). (2.5) 

Then, the wave profile r] is real- analytic r] G C"^(S). 

Remark 2.1. It is clear from the proof of Theorem 12.11 that its assertion is true in the case 
£7 = for arbitrary analytic vorticity functions 7,7, provided (I2.5p is fulfilled. When a > 0, 
and 7,7 G C"^(M) are bounded and satisfy additionally 

7' < and 7' < 0, (2.6) 

then the claim of Theorem 12.11 is still true. Relation (|2.6p ensures the unique solvability 
of the Dirichlet problems (|2.4ap . see Theorem 3.3 in |24| . which permits us to express the 
problem (|2.4p by the equation ()2.7p . Our method may be used to prove analyticity of the 
profile for traveling water waves with stagnation points in the capillary and capillary-gravity 
regime, provided the vorticity function is non-increasing, cf. |22|I23|. 

Proof. Assuming cr = first, we re- write the coupling equation (|2.4bp in a different way. 
Since the boundary conditions ony = rj[x) of both problems ()2.4ap imply that the tangential 
component of Vij: and V'0 vanish at the wave profile, this implies 

IVV'P = (S^V")^ and \Vi^\'^ = {dJ^f on y = r]{x). 

Here, 1/ := {—rj', l)/(l-|-?7'^)^/^ is the unit normal at y = r]{x). With this notation, equation 
()2.4bp is equivalent to the following relation 

f{x,y,d^il;,dui^) = p{d^tljf - p{d„tl;f + 2g{p - p)y -Q = on y = r/(x). 

Clearly, / is real-analytic in all its arguments. Furthermore, our assumption (|2.5p implies 
that 



ony = r]{x), which shows that the assumptions of Theorem 3.2 in |26| (see also the Remark 
following it) are all satisfied and the claim follows at once. 
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The proof in the case when o" 7^ is different. Indeed, we can reduce the problem ([2 
(see e.g. Section [3]) to an operator equation 

+HV) = 0, (2.7) 



(1 + 7/2)3/2 

where setting 

U:={r]e C2+°(S) : < 1}, (2.8) 

we have that (j) -.U -f C^+"(S) is a real-analytic operator with (/)(C°°(S)) C C°°(§). Because 
we use parabolic theory later on in the proof, we introduce the small Holder spaces /i'^+"(S), 
A; G N, as the closure of the smooth functions C°°(S) in C'^^"(§). We define now D := 
U n /i2+"(§), and observe that ^J}) may be written as 

^'=(17^+^^^)' ^^-'^ 

with ip G C'^{D, h^^"{S)) given by ip{r]) := r]' + (p{r]), t] £ D. Indeed, one can easily see from 
([2:91) that rj £ C°°(S), which imphes r] e D. 

Let now ^ : M ^ /i2+"(S) be the function defined by ^(i, x) = r]{t + x).l{ we remark that 
(p{i]{t + •)) = ip{i]){t + •) for all i £ M, which is a direct consequence of the unique solvability 
of (|2.4ap and of the fact that the variable x does not interfere into ()2.4ap , we find 

(1 + (^xO )^ 

:it, x) + ^{v{t + ■)){x) = —-^i-—{t,x) + mt)){x). 



{i + {d.CW/' ' (1 + (5.0^)^/2 

Whence, ^ is a solution of the autonomous problem 

dt^ = ^0, i>0, m=r], (2.10) 
where <I> G C"^(Z), /i°(S)) is the mapping 

Given ^ E we have a^>(e)[C] = (1 + {d^O'^)~'^^'^d^C + H with A G £(/i2+°(S), /ii+"(S)). 
Using well known interpolation properties of the small Holder spaces: 

(/i'^o(§),/i'^i(S))e = /i(i-^)'^"+^'^i(S) 

if G (0, 1) and (1 — 9)aQ + 9ai £ M"*" \ N, where (•, ■)g denotes the continuous interpolation 
method of DaPrato and Grisvard |12| ( see also [Il|28j ), we obtain from Corollary 3.1.9, 
Proposition 2.2.7, and Proposition 2.4.1 in |28) that the Frechet derivative <9'5(^) is the gen- 
erator of an analytic semigroup in i2(/i"(§)). Corollary 8.4.6 in |28| ensures that the unique 
solution of ([230]) to the initial data ^(0) = 77 is analytic, that is ^ G C"^((0, 00), /i2+"(§)), 
which implies the desired assertion. □ 
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Internal water waves with a critical layer 7 

3. Bifurcation analysis and the main result 

In the remainder of this paper we restrict our analysis to the stable regime when the fluid 
in is more dense than that located above, that is p > J), and the fluids have constant 
vorticity 

a; = 7GM and cJ = 7 e M. (3.1) 

Furthermore, the surface tension coefficient may take any value a € [0, 00). 

As a first step, we introduce a constant A into the problem which will allow us to describe 
the laminar fiow solutions of ()2.4p . i.e. sohitions with a fiat wave profile, located at y = 0. 
Later on, we use this constant as a bifurcation parameter to find non-fiat solutions of ()2.4p . 
To this end, when r/ = 0, we observe that the functions (i/j, tp) =: [ipQ, ^q) solving ()2.4ap are 
given by 

2 2 

ipo{x,y) ■■= ^ in and i^o{x,y) := ^ + Ay in Q^, (3.2) 

provided the constants m and m satisfy 

m = uj/2 and m = A+ZI;/2. (3.3) 

Hence, given A € M, the tuple [r] , ijj , il^) := (0, "00) V^o) ^ solutions of (|2.4p if the constants 
m,rn, and Q are given by (|3.3p and (|2.4cp . 

In order to determine non-fiat solutions of ()2.4p we re- write problem ()2.4p as a nonlinear 
and nonlocal operator equation having (A, rj) as unknowns. Because the domains where the 
Dirichlet problems (|2.4ap are posed depend upon r], we need to transform these problems 
and the equation ()2.4bp on fixed reference manifolds. Therefore, we define the functions 
(/)^ : — 7- Q,^^ and 0^ : — )■ fi* by the relations 

(pn{x,y) := {x,y+{l + y)r]{x)) hi and ^^(x, y) := (x, y (1 - y)r/(x)) in f^f,, 

and observe that if rj belongs to U (see relation (|2.8p ). then (j)^^ and (j)^ are diffeomorphisms. 
We use these diffeomorphisms to transform the Laplace operator into a differential operator 
on the fixed domains fig and fig, respectively. More precisely, setting 

Airi)w:= {Aiwo^-')) o^^, wgC^+''{%), 
we find the following expressions for the differential operators A{r]) and A{rj) : 
AM-d^ -2^1+^9^ + ' -(l + y)il±^^>tV!5 

•^Kv)-Cxx ^ cixy^y (i-??)2 {i-rjyj ^ ' (i-??)2 



for sl\ rj ^ U. >From these explicit expressions we can easily see that the functions r/ 1— t- 
A{i]) and 1] I—)- A{ri) are both real- analytic. Furthermore, corresponding to the coupling 
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condition (|2.4bp we define the boundary operators B : U x C^+"(r2Q) — > C^+°(S) and 
B:Ux C2+"(0^) C^+°(S) by the relations 

S(7?,zz;) :=tr(|V(u;o^-i)|2o,/,^) and ^(r/, u;) := tr (|V(lZJ o o , 

with tr being the trace operator with respect to the hne S x {0} = S. For these operators 
we find 

B{r,, w) := tY{d^wf - tr d^^w tr dyW + 7^^^ tT{dyw)\ (77, w) (^U x C^+^ , 

1 + 7/ (1 + r/j^ 

B{r^,w) := ti{d,w)^ - ,^tvd,wtvdyW + ^^^!^tT{dywf, {7],w) x C2+°(0^), 

which shows that B and B also depend analytically on their variables. 

To conclude, we observe that if ■0 and ij) are the solutions of the problems (j2.4ap (for 
some r] £ U) when lj,uJ are constant and the constants m,rn are given by (j3.3p . then 
HJ := tp o (pjj =: T{'q) and W := ip o (j)^ =: T{X,r]) are the solutions the solutions of the 
Dirichlet problems 

' A{r])w = Lo in Qq, ( A{r])w =lJ in Qq, 

w = on y = 0, and < w = on y = 0, (3.4) 

w = LO /2 on y = —1 ItJ = A + aJ/2 on y = 1, 



respectively. Since, A and A are real-analytic, and the right-hand side of the equations (|3.4p 
depends analytically on A, we obtain that T and T are also real-analytic. For the proof we 
refer to Lemmas 2.2 and 2.3. in |17| . 

With this notation, finding the solutions of the coupled problem (j2.4p when 

the fluids have constant vorticities and the constants m^rn^Q are given by relations (j3.3p . 
(j2.4cp . reduces to determining the solutions (A, ?y) of the nonlocal and nonlinear equation 

^(A,r?) = 0, (3.5) 

where ^ -.^xU ^ C"^'+"(S) is given by 

^'(A, 77) :=pfi(r/, T(A, 7?)) - pB(7?, T{ii)) + 2g{-p - p)r^ + 2ar?7(l + 7?'2)3/2 

{pB{v, T(A, 7?)) - pB{7], T{v)) + 2g{p - p)rj{x)) dx, (A, r/) G M x ZY, 



with cr' = if (7 > and a' = 1 for ex = 0. 

We note that the laminar flow solutions found at the beginning of the section are in 
correspondence with the trivial solutions (A, 77) = (A, 0) of ()3.5p . By applying the bifurcation 
theorem from simple eigenvalues due to Crandall and Rabinowitz |11) to equation ()3.5p we 
show in Theorem 13.11 that infinitely many analytic branches consisting of non-flat solutions 
of (|2.4p intersect this trivial set of solutions of (|3.5p . To this end, we restrict first the 
domain and range of This is done by introducing the subspaces C'l^Q^yi'S') , k,m G N, 
of C"^"'""(S) which contain only even, 27r/k periodic functions with integral mean zero. 
Choosing functions with integral mean zero corresponds to a choice of equal volumes of 
both liquids beneath the rigid lid, meaning that, at rest, the flat wave profile is located at 
y = 0. We use next the intrinsic properties of (|2.4p to prove the following lemma. 
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Internal water waves with a critical layer 9 
Lemma 3.1. Let a G [0,oo) andlAkfi^ev •= ^ operator'^ is real-analyti 

M/GC"^(Mx^/fc,o,e.,<ote°(S)), 

and 

2ar]" 



^{\, 1]) =pB{rj, T{X, v)) - pB{7], T{r,)) + 2g{-p - p)r^ + 



(1 + r/'2)3/2 
for all {X,T]) G M X Ukfi,ev 

Proof. The analyticity of ^ is a consequence of the fact that the operators T,T,B,B are 
real-analytic in their variables. 

Let now (A, ry) G M x Ukfi,ev be given. Since t] is even, we may define the functions 
wi{x,y) := w{—x,y) for {x,y) G and wi{x,y) := w{—x,y) for G Qq, where w and 

It; denote the solutions of ()3.4p . respectively. Since the Dirichlet conditions are constant, it 
may be easily checked that wi and wi are solutions of the problems ()3.4p . respectively, so 
that, by the weak elliptic maximum principle, w = wi and W = Wi. A similar argument 
shows that w and W are both 271 /k periodic in x and we obtain, by using the explicit relations 
for the boundary operators B and B, that ^{X,rj) is 27r/k periodic and even. To finish the 
proof, we note, with our choice of Q in (|2.4cp . that ^{X,r]) has integral mean zero if rj has 
this property. □ 

An important step in our analysis is to determine the Frechet derivative of ^ at the 
trivial branch of laminar solutions. 

Lemma 3.2. Given A G M, the Frechet derivative dr,^{X,0) G >C(C^+|'g^(S), C^q+'^(§)) is a 
Fourier multiplier 

oo oo 

(9^^(A, 0) Ofc cos(A;x) = fik{X)ak cos{kx) (3-6) 

k=l k=l 

with symbol 

k 



g{p- p) -ak'^ +ujpX + P ^^^^j^^ 



k>l. (3.7) 



Proof. Taking into account that ^(O) = ipQ and ^(A, 0) = ipQ for all A G M, we obtain from 
9^i3(0,Vo) = 0, d^B{0,i;o) = 0, 



dr,B{0, V^o) N = 2A2r7, d^B{0, iP^M = 2A tr dyW, ^^'^^ 
and the chain rule that 



drj^{x,o)[i]] =2px ^tT{dyd^T{x,oM) - j^ti{dyd^T{x,oM)dxj (3.9) 

+ [2g{p - p) + 2X^p]r] + 2ar]" (3.10) 

for all 7] G C^+'^^XS)- We determine now the Frechet derivative w := drj'T{X,0)[rj]. Therefore, 
we differentiate the equations of the Dirichlet problem (|3.4p solved by 7~(A, rj) with respect 
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to rj at 7] = 0, and find that IZJ is the solution of the problem 

Aw = -dA{0)[r]l^o in fl^, w = on dn^, (3.11) 

where, by the explicit expressions found at the beginning of the section 

dJiOMlPo = 2oJr, - (1 - y){u]y + A)r/". 

In order to determine relation (|3.6p we use Fourier expansions for functions in q (^) 
make a similar ansatz for the solution w := driT{X, 0)['i]] of (|3.1ip : 

oo oo 
rj = Ofc cos(A;x) and w = akuikiu) cos{kx). 
k=l k=l 

The right-hand side of the first equation of (|3.1ip may be then expanded as follows 

oo 

-dA{0)[r]]ipo+ = akbkiy) cos{kx) with bk{y) := -2lJ - (1 - y){uJy + X)k'^, 
k=l 

and, plugging all these expansions into (|3.1ip and matching the coefficients corresponding 
to cos(A:x), yields that Wk solves 

w'l - = bk, < 2/ < 1, and Wk{0) = Wk{l) = 0. (3.12) 

The solution of (|3.12p is given by 

Wk{y) = ^^^^^^^ sinh(A;y) - A cosh(A;y) + Lj{y - y^) + A(l - y), (3.13) 

and together with (j3.9p we obtain the desired relations p.6p and p.7p . □ 

We recall now a theorem which provides a sufficient condition for an operator in order 
to be a Fourier multipliers between Holder spaces of periodic functions. 

Lemma 3.3 ( |18l Theorem 3.4] ). Let r,s be two positive non-integer constants and let 
{Mp)p^z C C be a sequence satisfying the following conditions 

(i) sup |p|''"'^|Mp| < oo, 

p6Z\{0} 

(ii) sup |pr-"+i|Mp+i - Mp\ < oo, 
pez\{o} 

(m) sup \p\''-'+'^\Mp+2 - 2Mp+i + Mp\ < oo. 

pGZ\{0} 

Then, the Fourier multiplier satisfies 

^Opi'P^ ^^Mpapi'P^ G £(C^(S),C"(S)) 
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Internal water waves with a critical layer 11 

Using Lemma [3.3l and the relations ()3.6p and (j3.7p . it is not difficult to see that drj^{X, 0) 
is an isomorphism if A is chosen such that /ifc(A) 7^ for all integers > 1. To find these 
values of A we solve the quadratic equation //^(A) = and find the solutions 



u} tanh(A;) , / q(p — p) + ak"^ tanh(fc) u"^ tanh^(A;) 

^-=-2^+(-^)V^^^i fe^ + T^;^' ^ = ^'^- 

We consider the case a = first. Since x 1— t- tanh(x)/x is a strictly decreasing function 
mapping [0, 00) onto (0,1], and, for > 0, t 1— )• t + + Bt is increasing (resp. t 1— )• 
t — + Bt is decreasing) on [0, 1], we conclude that (A^) are strictly monotone sequences. 
If cr 7^ 0, we choose do > such that (A^) are strictly monotone sequences for all a > a^. By 
Lemma 13.31 we see that in both cases (?r;^(A^, 0), i = 1,2, is a Fredholm operator of index 
zero having a one dimensional kernel. 

This leads us to the existence result of this paper. Besides existence of analytic curves 
consisting of traveling internal waves, we determine the second order Taylor expansions for 
these curves in a neighbourhood of the laminar flow solutions. These expansions give us 
sufficient information to provide the precise picture of the streamlines (the level curves of 
the stream functions) in the frame moving with wave speed c, see Theorem 14.11 

Theorem 3.1. Assume that < p and if a ^ 0, then let o" > (Tq. Given A; G N, /c > 1, there 
exists Ek > and real-analytic curves 

{XI, Vl) : {-ek,ek) ^ M x C2+"^,(§), i = 1,2, 

consisting only of real- analytic solutions of (j3.5p of minimal period 2'ir/k and having exactly 
one crest and trough per period. These are the only solutions of (|3.5p close to (A^.,0), and 
for s — 7- UJe have 

Xl{s) = K + 0{s% ^^^^ 
r]i{s) = -cos{kx)s + aicos{2kx)s'^ + 0{s^) m C^+"^JS), 

with constants a\ given by (|3.23p (with A replaced by A^J. 

Proof. We verify flrst that the assumptions of the theorem on bifurcations from simple 
eigenvalues due to Crandall and Rabinowitz |1H Theorem 1.7] are satisfied. We already 
know that, when A = A^ for some /c G N and i= 1,2, the Frechet 9^>I'(A^,0) is a Fredholm 
operator with 

Ker a;,^'(A*fc, 0) = span{cos(A;x)} and Im 5^1- (A^, 0) span{cos(/cx)} = C^'o^e" (§)• 
Furthermore, diff^erentiating (|3.6p with respect to A we obtain that 



I g{p — p) -\- ak"^ k u? 

tanh(A;) 4 



9|^^(Al,0)[cos(fcx)] == ±4pJ ^^^ py^^^ +^cos{kx) 



which implies 5|^^'(A^, 0)[cos(/ca;)] ^ Im(?^^'(A^, 0). The existence of the analytic bifurca- 
tion curves follows now from the above mentioned theorem and Lemma [3.11 

We pick now a solution {X{s),r]{s)) of (j3.5p located on one of the curves (A^, ry^), i = 1,2, 
and denote by tj; := T{r]{s))o(j)'^^^^ and ■0 := ^(A, 'n{^))°4'ri{s) stream functions associated 
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12 A.-V. Matioc 

to it. In order to prove that ri(s) is real-analytic, we show that the assumption ()2.5p is fulfilled 
provided e/j is sufficiently small. Indeed, since ^(O) = "00 a-nd T(A^, 0) = ipQ, with -00 and ip^ 
given by (|3.2p . we obtain from jVV'oP + I^V^oP = |Ay^>Oony = that (|2.5|) is satisfied 
by the laminar flows (A^(0), 7/^(0)), i = 1,2. Choosing small enough, the real-analyticity 
of r](s) follows now from Theorem 12.11 by making use of the continuity of the bifurcation 
curves and of the solutions operators T and 7~. 

Further on, we prove the asymptotic expansions ()3.14p and show that the internal travel- 
ing waves we obtain have exactly one crest and trough per period. To this end, we fix E N, 
A; > 1, and, to ease notation, we let in this final part of the proof A := A^,, for i £ {1,2}, 
and denote by (A,7/) the corresponding branch of solutions (A^,ry^). Since the bifurcation 
curve (A,?/) is real-analytic, we obtain from Theorem 1.7 in |11| that 

A(s) = A + A,(0)s + O(s2) and ij{s) = -s cos{kx) + t{s) + 0{s^) mCl+1^{§), 

(3.15) 

with r(0) = Ts{0) = and where the index s denotes the derivative with respect to 
the variable s. Additionally, the function r takes values in the closed complement Xq of 
span{cos(A;x)} in Cfc^e„(S)- Proceeding similarly as in |34) . we find first from (|3.15p that 

r]' = sksm{kx) + 0{s'^) and r]" = sk"^ cos (kx) + 0{s'^) in C(S). 

Therefore, if'iO) > 0, ?/'(7r//c) < 0, and, since ?/ is odd, we also have r?'(0) = r]'{TT/k) = 0. 
We resume that ij' is positive on (0, vr/A;) provided Sk is small, meaning that the wave has 
its creast located at x = n/k and the trough at x = 0. 

The next step of the proof is to determine the derivatives A^ (0) and Tgs (0) . Differentiating 
the relation ^'(A(s), t/(s)) = twice with respect to s we find, at s = 0, that As(0) and Tss{0) 
are related by 

- 2A, {0)dl^^{A, 0) [cos(A;x)] + d^^^{A, 0) [cos{kx)f + d^^{A, 0) [r,, (0)] = 0. (3.16) 
Clearly, we need to find the second order derivative 

d^^'^{X,0)[cos{kx)f =p (l - j 7dx^ -p(^- j I dx^ , (3.17) 

where, setting ^ := cos(A;x), we made the following notation 



/ :=a,\^(0, Vo)[^]' - 2d'^^BiO,i;omdTiom] - ^^^(0, Vo)[ar(0)[e]]2 
- ^^BiO,^Po)[^'T{0)m, _ 

I :=d'r,r,m + 2d^,^m %m d^nA, o) m + ^^^(o, ?o) [5,r(A, o) 
+d^^o,i^o)[dLnA,om% 



(3.18) 



In view of (|3.2p . we compute 

d^^m^o)[e = 2A2 (3^2 + ^.2) ^ d^^^BiO,^Q)[C,v] = AACtrdyV - 2i'Ai,d,v, 

dl^B{ri,w)[vf = 2tr {{djvf + {dyvf) , dl^B{r^,w)[vf = 2tr {{d^vf + {dyvf) , 

a2^B(o,Vo) = 9^^mV'o) = o, 

and, recalling (|3.8p . we have (9u,;B(0, ■(/'o) = 0. Consequently, we need to determine only the 
derivatives <9T(0)[^] and d'^^T{A,Q)[^f to obt ain an explicit expression for the right-hand 
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Internal water waves with a critical layer 13 

side of (j3.17p . Concerning dT{0)[(,], we have to study a linear Dirichlet problem similar to 
p.lip . and one finds 

dr{0)[cos{kx)] = a;(y2 + y) cos(fcx). (3.19) 

As for the the second order derivative, we differentiate the Dirichlet problem (j3.4p for w 
twice with respect to r/ and obtain, at r/ = 0, that v := d'^^T{A, 0)[^]^ is the solution of the 
problem 

Av = -2dA{0)[^]dr^T{A,0)[^]-d^A{0)[^,^]i^o mQ^, v = ondnl, (3.20) 
where 

d^Amie^ = uj (2(1 - yfe + Qe) - (i - y)i^y + a) (4^" + 2^e) , 

dA{0md^T{A, om = (-2(1 - y)^'dly + 2idly - (1 - y)edy) {d,T{A, 0)[e]) . 

By (j3.13p . we have 9^T(A, O)[cos(fcx)] = Wk{y) cos{kx), so that we can express the right- 
hand side of the first equation of (j3.20p as follows 

2dA{0md^T{A, om + d^A{om^i^o = Eo{y) + E^kiy) cos{2kx), 

whereby 

EM - + 2t^A - CO*,..)) - (1 - .).»A - .nK.„ 

/ sinh{ky] 
\ tanh(A;) 



E2k{y) := - cJ + 2A;2zj(l - yf + 2k^A ( ^^^^^^^ - cosh(%) 



+ 3(1 - y)k^A - sinHky) 



Due to the linearity of (|3.20p . we write d^^T{A, 0)[cos(A;x)]^ = Ck{y) + dk{y) cos{2kx), with 
Cfc denoting the solution of ()3.12p when (k, bk) = (0, —£"0) and being the solution of ()3.12p 
with (k, bk) are replaced by {2k, —E2k)- Since all Frechet derivatives of 7~ and T with respect 
to T] have zero boundary values, of relevance for our purpose are only the first derivatives 



4(0) :=-fe'A-| 

,/ /^x ,2a ^A 2k'^A _ kZo 

4(0) := A — — H — — TTTTTT + 



tanh(/c) taiih{k) tanh(2/c) tanh(2A;) 
Summarising, we have shown that 

/ =uj'^ + uj"^ cos{2kx), 



-r f kA Y 2kA^ 2 a a2,2 



\tanh(fe) J tanh(A;) 



+ 



' +u}] + — — — + — — - 3k^A^ 



ytanh(A;) J tanh(A;) tanh(2A;) tanh(2A;) 



cos(2/cx). 
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14 A.-V. Matioc 

Recalling (|3.17p . we find that d^j^^{A,0)[cos{kx)]'^ = AkCos{2kx), whereby 
' / kA V 4A;2A2 2fea;A 



A, 



\tanh(/c) 



+ uJ + 



+ 



tanh(A;) tanh(2A;) tanh(2/c) 



(3.21) 



To finish, we observe that (9^^^'(A, 0)[cos(A;x)]^ belongs to Im drj"^ {A, 0) , meaning that 
cos(A;x) is orthogonal on (9^^^'(A, 0)[cos(A;x)]^ in L2(S)- Whence, if we multiply the relation 
(j3.16p by cos{kx) and integrate it then over the unit circle, we get that Xs{0) = 0. Moreover, 
since the restriction dri^{A,0) : Xq — )• Im9^^(A, 0) is an isomorphism and Tss{0) £ Xq, 
(TOD yields 



TssiO) 



-{d^^{A,0))-^d^^^{A,0)[cos{kx)]\ 



Setting 



Oik ■-- 



2/X2fc(A)' 



(3.22) 



(3.23) 



we conclude from (|3.6p that Tss{0) = 2akCOs{2kx), and together with (|3.15p we find the 
desired expansion for i]. This completes the proof. □ 



4. Streamlines for internal waves with stagnation points on the profile 

In the remainder of this paper we restrict our considerations to the curve (A^, ry^) and denote 
by (A, rj) a solution of (j2.4p which lies on this curve. As in the proof of Theorem 13.11 we 
write -0 and ip for the stream functions associated to this solution. The next theorem shows 
that the traveling wave solutions found in Theorem 13.11 possess points which are stagnation 
points when considering the wave profile as a part of the fluid located below, but loose this 
feature when considering the interface as being a part of the fluid located above. Particularly, 
as we approach these points, the velocity (u, v) of the fluid particles located below the wave 
profile satisfies {u,v) — )• (c, 0). 

We shall exemplify this in Theorem 14.11 in the particular case when tJ > and a; > 0. 
We may choose also u < 0, and only the orientation along the streamlines beneath the wave 
profile has to be changed in Figure [T] Allowing uJ to be negative, we may obtain stagnation 
points within the fiuid located above. However, the picture of the additional critical layer 
we could obtain in ri* has been studied in detail in [91I191IM]. 

Theorem 4.1 (The streamlines in the moving frame). Additionally to Theorem \3.1[ 
assume that cJ > and a; > 0. Then, within a period, the streamlines corresponding to a 
solution {X,ri) on the curve (A^.,?7^.) are qualitatively described in FigureUl Particularly, there 
exists a critical layer consisting of closed streamlines which is delimited from above by the 
wave profile and from below by a separatrix which connects two stagnation points (solutions 
ofVip = 0) which are located on the wave profile. Furthermore, there exists exactly one more 
stagnation point which is situated in the center of the critical layer. 



In order to prove Theorem 14. H we rely on the expansions (j3.14p . Compared to 
the situation we consider is more degenerate, because, by our choice (|3.3p of the constants, 
bifurcation occurs when ipQ^y = on y = 0, and we need second order expansions for the 
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Fig. 1. This figure illustrates the streamlines in the moving frame in the case when > 0, U > 0, and 
A = A^. The dashed blue curve is the critical level where i/iy vanishes and the dashed red curves consist 
only of points where vanishes. The critical layer of closed streamlines is located between the wave profile 
and the separatrix (the thick line below the wave surface) which connects the two stagnation points on the 
wave profile. A further stagnation point is located below the wave profile on the line x — n/k. 

bifurcation curves, cf. ()3.14p . to be able to analyse the streamlines within the fluid domains. 
First, we prove: 

Proposition 4.2. Let (A, ?y) be located on the bifurcation curve (A^,??^), and define 

Ofo := {{x,y) G : < x < vr/A;}. 
If Ek is small and to ^ 0, we have the following properties: 

(a) IfuJ> 0, then -i/'^ > in fi^^ and ip^^x, y) < for all x G (0, n/k) and rj{x) < y < 1; 
(6) ujil^yy > in 0,^; 

(c) There exists a unique point ( G (0,7r/A;) with C = vr/(2A;) + 0{s), and a continuous 
curve y(^ : [0, 27r/A;] — )• M which satisfies 

(i) yc_{x) = r]{x) for x ^ (C, 2-K/k — C) and —1 < yQ{x) < r}{x) for all x G (C, 271 /k — 

C); 

[a) ujipy > for X G (C,27r//c — () and y G {y^{x),r]{x)], ipy{x,yi^{x)) = for all 
X G [C, 27r//c — <^], and ujipy < elsewhere in Q^; 

(d) There exists a curve x = ^(y) : [—l,yb] (0,7r/A;) such that iS,iyb),yb) is o, point on 
the wave profile and 

' ujipxy{x,y) > 0, ify < yb and x < ^{y); 

< uj'ipxy{x,y) = 0, if X = C{y); (4.1) 
^ L0'i/Jxy{x,y) < 0, ifyb <yor[y<yb and ^(y) < x < ir/k]. 



Moreover, ^ = 7r/(2/c) + 0{s) uniformly in y & [—l,yb]. 
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16 A.-V. Matioc 

(e) There exists a curve x = ^(y) : [— l,y{,] — )■ (0, vr/A:) such that ('^(yf,), ^f,) is a point on 
the wave profile and 

LOTpxix^y) >0,ify<yp and x < ^(y); 

ujip:,{x,y) = 0,ifx = C{y)] _ (4.2) 

ujtpx{x,y) <0,ify>yh or [y <yf, and i{y) <x < Ti/k]. 

Moreover, (C, ^((C) = ^h) and^ = tt / {2k) + 0{s) uniformly zn y G [-1, y^,]. 

Proof. To prove the claim (a), we recall that ijj^^y = 2iJy + > in fig. Using the 
continuity of the operator T and of the bifurcation curve (A^,??^) it follows that > in 

O*. Furthermore, we notice that ip^{x,r]{x)) = —ri'{x)iljy{x,rj{x)) < for all x £ (0,7r/A;), 
and since ■0 is even in x we find that -0^ < on the boundary d^t, whereby ilf '■= {{x, y) € 
0* : X G (0,7r/A:)}. Recalling ()2.4ap 9. we find that Aip^ = in and the claim (a) follows 
from the strong elliptic maximum principle. 

We consider in the remainder of the proof only the fluid located below the interface r]. 
The functions rj and ip being both even functions, it suffices to restrict our considerations 
to the domain Qj,. In order to prove (b) and (c), we use the first order Taylor expansion 

2 

w := r(r?) = r(0) + dTiOM + 0(s2) = ^ - sdTiO)[cos{kx)] + 0{s^) in C^+^'i^). 



Recalling ()3.19p . we get 



w 



ujy 



2 



suj{y'^ + y)cos{kx) + 0{s'^) in C72+"(0^). (4.3) 



Since ip{x,y) = w{x,{y — 7]{x))/{l + ri{x))) for all (x,y) G we immediately obtain 
the assertion (6). Furthermore, at y = —1, we have that dyw{-,—l) = —uj + 0(s), while 
differentiating ()4.3p at y = 0, yields 



dyw{x,Q) = —sujcos{kx) + 0{s'^). 

Repeating the arguments in the proof of Theorem 1 3. 11 we find that ujdyW is strictly increasing 
on [0, 7r/fe]. Since u}dyw{0, 0) < and ujdyw{Tr/k, 0) > 0, we deduce that there exists a unique 
C G (Ojir/k) such that dyw{C,0) = 0. By virtue of (6), we find for each x G (Cvr/A;] a point 
y^(x) G (—1,0) with the property that ujdyw{x,y) < for y < y^(x) and ujdyw{x,y) > 
for y > y^{x). The mean value theorem shows that is close to 7r/2k in the sense that 
C = TT/{2k) + 0{s). We appeal now to the fact that dyW is even to obtain the desired 
conclusion (c). 

In order to prove (d) and (e), we are obliged, to use a better approximation for T{ri) 
w =T{r]) = r(0) + dT{OM + ^d^T{0)[vf + 0(8^) 



^ - sdT iO)[cos{kx)] + (dT{0)[ak cos{2kx)] + ^d'^T{0)[cos{kx)Y j + 0(s% 

_ (4.4) 
in C2+«(0[;), cf. dSTlD. To determine d^T{0)[cos{kx)]'^ we differentiate the Dirichlet problem 
corresponding to fij^ twice with respect to r/, and find that d'^T{0)[cos{kx)]^ is the 
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solution of the Dirichlet problem 

Av = -25^(0) [cos(x)]ar(0)[cos(a;)] - d^A{0)[cos{kx)fipo ^^^o, v = ondQ^o- 
Similarly as in the proof of Theorem I3.H we find that 

d'^T{0)[cos{kx)f = '^^y'^^y^ + ujh{y) cos{2kx), (4.5) 



whereby 



For the clarity of the exposition we leave the proof of (|4.5p to the interested reader. 
We sum now the relations (|4.4p - (|4.6p and use (|3.19p to conclude that 



LO 2 



+ afe(/ + y) + ^^ cos(2A:x 



(4.7) 



in C'^^°'{Qq). The next goal is to find the expansion corresponding to ipx- Recalling (|3.14p . 
we have 

^ r]'{l - 2??) + 0{s^) = sk sm{kx) + s^k{l - 2ak) sm{2kx) + 0{s^) 



(1 + r?) 

= y{l-r])-7] + 0{s^) = y + {y + l)scosikx) + 0{s^), 

l + r] 

and taking into account that ip{x, y) = w{x, {y — r])/{l + r])) we find from 
ipx{x,y) Wx ( y-r]{x)\ {I + y)r]' Wy ( y - 7]{x) 

— X, ; -t; X 



u) (x! \ ' 1 + r/(x) J (1 + 7^)2 (x! \ ' 1 + r/(x 

that 

=^ (SS!m ^ + 0{s') in C2+"(^). (4.8) 

This is the key point in the proof of (d) and (e). 

To keep the notation short we introduce the auxiliary function / : [— 1, 0] — )■ M defined 

by 

, A; /sinh(2A;y) , A sinh(2/c(l + y)) , ^ 

Since tpxy = s'^ujf'f^{y) sin{2kx) -\- 0{s^) and is strictly positive in [—1,0], we immedi- 
ately obtain the assertion (d). Additionally, for arbitrary 6 G (0, -7r/(8A;)), the derivative 
^''Pxxyi^^y) < for all {x,y) G with x G [7r/(4A:) + (5, 37r/(4/c) — 5], provided is suf- 
ficiently small. Recalling that ■0 = uj/2 on y = —1, we get tp^xix, —1) = and therefore 
^'4'xx{x,y) < for all x G [7r/(4A;) + (5, 37r/(4fc) — 5] and y > —1. On the other hand, the 
mixed derivative coipxyi'^ / i^k) + 5,y) > and a;'(/'a;j/(3vr/(4A:) — 6,y) < for all y, and taking 
into account that ^^^(x,— 1) = 0, we conclude that there exists a point G (minr/, maxr]) 



February 15, 2012 1:27 WSPC/INSTRUCTION FILE JNMP_Matioc_Anca 



18 A.-V. Matioc 

and a curve x = ^(y) : [— l,y{,] — )■ {0,Tr/k) such that ('^(yf,), yf,) is a point on the wave profile 
and 

V'x(?(y),y) =0, for ahye [-l,y,]. (4.9) 

In fact Ub = yg(0, and ^(y) = TT/{2k) + 0{s) uniformly in y G [— l,yj]. This curve [x = ^(y)] 
splits the domain $7;, into two subdomains \ [x = (,{y)] = Q^jU 51^, being the domain 
which has x = as boundary component. Since Ail^x = in Oj^ and, since by (c), ujipx > on 
5r2[,, the strong elliptic maximum principle ensures that uipx > in The same argument 
shows that utpx < in ^2^. 

Finally, in order to show that {Cy^iC) = ('^(y6)>yfe)) differentiate the relation 
ip{x,'q{x)) = with respect to x and obtain, by virtue of r;' > in (0,ir/k), that, if 
X € (0, tt/Zc), then ipxix, y) = if and only if ipyix, y) = 0. We infer from (c) that the desired 
equality holds and the proof is completed. □ 

We come now to the proof of Theorem 14.11 It is based on Proposition 14.21 and the fact 
that the curves obtained in Proposition 14.21 (c) and (d) never intersect. 

Proof of Theorem 14. IL In order to precisely determine the streamlines within the fluid 
located below we need to relate the two points (Ciyc(C) and (■^(yb), yt). Therefore, we dif- 
ferentiate the equation ip{x,ri{x)) = twice with respect to x and find, at x = that 



V 1 



CO 



This relation is obtained by also using the relation ^y(C,y^(^)) = together with the 
equation A?/^ = a; in Q^, cf. (I2:iall . Moreover, by (f3J5]l and (OD . we find the following 
expansions 

r,'2 ^ ^2^-^ _ cos(2A:2;))s^ + ©(s^), (4.10) 

^ = k'( + cosh(2A:y)^ cos{2kx)s^ + 0{s'), (4.11) 
CO \ tanh(2K) J 

which yield, in the end. 



id 



/ sinh(2fcy) 

V tanh(2/c) ^ 



cos(2A;C) ^s^ + 0(s^). 



Taking into account that (C5yc(C)) ^ ('''"/(^fc); 0) and r/'(C) > 0, we conclude that 
and, by (liTll . 

Vb < yc(C)- (4.12) 

This implies that the curve [y = y^] is located entirely in the region [oJil^xy < 0] and that y^ 
is strictly increasing. Indeed, from the relation 'ipy{x,y(^{x)) = we find that 

'fpxyix, y({x)) + y'(-{x)^yy{x, yc_{x)) = 

for all X G [C,27r//c]. Since i^xy{C^yc,{C)) < and tpyy > 0, we conclude that y^ is a strictly 
increasing function on [^,27r/A;]. Whence, for x G (Cn/k), the point (2;,y^(x)) is located in 
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the region where ipx < 0, cf. (j4.2p . This means that the function [?], vr/A;] 9 x i— t- uj-ip{x, y({x)) 
attains its minimum in x = ir/k and, by (c), this value is also the minimum of ojip : 

mm^J^l^ = u)tlj{Tr/k,y(^{'ir/k)). 

We resume that there exist exactly three stagnation points one period: two on the wave 
surface (Cj 2/^(0) c^iid (27r//c — C,y^{27r/k — C)), and one (vr/fc, y^(7r/A;) beneath the wave 
crest. Since by Proposition 14.21 (c) and (e) we know the sign of ■z/'x and ipy in the whole 
domain $7^ , we conclude that the streamlines of the flow are as in Figure [TJ □ 



References 

[1] H. Amann: Linear and Quasilinear Parabolic Problems, Volume I Abstract Linear Theory, 
Birkhauser, Basel, 1995. 

[2] J. L. Bona, D. Lannes & J.-C. Saut : Asymptotic models for internal waves, J. Math. Pures 

Appl. 89 (2008), 538-566. 
[3] A. Constantin: The trajectories of particles in Stokes waves, Invent. Math. 166 (2006), 523-535. 
[4] A. Constantin and J. Escher: Analyticity of periodic traveling free surface water waves with 

vorticity Ann. of Math. 173 (2011), 559-568. 
[5] A. Constantin and R. S. Johnson: Propagation of very long water waves, with vorticity, over 

variable depth, with applications to tsunamis, Fluid Dynam. Res. 40 (2008), 175-211. 
[6] A. Constantin and W. Strauss: Exact steady periodic water waves with vorticity. Comm. Pure 

Appl. Math. 57(4) (2004), 481-527. 
[7] A. Constantin and W. Strauss: Pressure beneath a Stokes wave. Comm. Pure Appl. Math. 

63(4) (2010), 533-557. 

[8] A. Constantin and W. Strauss: Periodic traveling gravity water waves with discontinuous vor- 
ticity, Arch. Ration. Mech. Anal. 202 (1) (2011), 133-175. 

[9] A. Constantin and E. Varvaruca: Steady periodic water waves with constant vorticity: Regu- 
larity and local bifurcation. Arch. Rational Mech. Anal. 199 (2011), 33-67. 
[10] A. Constantin and G. Villari: Particle trajectories in linear water waves, J. Math. Fluid Mech. 
10 (1) (2008), 1-18. 

[11] M. G. Crandall and P. H. Rabinowitz: Bifurcation from simple eigenvalues, J. Fund. Anal. 8 
(1971), 321-340. 

[12[ G. DaPrato and P. Grisvard: Equations d'evolution abstraites nonlineaires de type parabolique, 

Ann. Mat. Pura Appl. 4 (120) (1979), 329-396. 
[13[ M. Ehrnstrom, J. Escher, and G. Villari: Steady water waves with multiple critical layers: 

Interior dynamics, to appear in J. Math. Fluid Mech. 
[14[ M. Ehrnstrom, J. Escher, and E. Wahlen: Steady water waves with multiple critical layers, 

SIAM J. Math. Anal. 43 (2011), 1436-1456. 
[15[ M. Ehrnstrom and G. Villari: Linear water waves with vorticity: Rotational features and particle 

paths, J. Differential Equations 244 (2008), 1888-1909. 
[16[ M. Ehrnstrom and E. Wahlen: On steady water waves with critical layers, preprint. 
[17[ J. Escher and G. Simonett: Classical solutions for Hele-Shaw models with surface tension. Adv. 

Differential Equations 2 (1997) 619-642. 
[18[ J. Escher and B.-V. Matioc: A moving boundary problem for periodic Stokesian Hele-Shaw 

flows. Interfaces Free Bound. 11 (2009), 119-137. 
[19[ J. Escher, A.-V. Matioc and B.-V. Matioc : On stratified steady periodic water waves with 

linear density distribution and stagnation points, J. Differential Equations 251 (10) (2011), 

2932-2949. 

[20[ K.R. Helfrich and W.K. Melville: Long nonlinear internal waves. Annual Review of Fluid 
Mechanics 38 (2006), 395-425. 



February 15, 2012 1:27 WSPC/INSTRUCTION FILE JNMP_Matioc_Anca 



20 A.-V. Matioc 

[21] D. Henry: Analyticity of streamlines for periodic traveling free surface capillary-gravity water 
waves with vorticity, SIAM J. Math. Anal. 42 (6), (2010), 3103-3111. 

[22] D. Henry: Analyticity of the free surface for periodic traveling capillary-gravity water waves 
with vorticity, J. Math. Fluid Mech. (2011), DOI 10.1007/s00021-011-0056-z. 

]23] D. Henry: Regularity for steady periodic capillary water waves with vorticity, to appear in 
Philos. Trans. R. Soc. Lond. Ser. A. 

]24] D. Henry and B.-V. Matioc: On the existence of steady periodic capillary-gravity stratified 
water waves, to appear in Ann. Sc. Norm. Super. Pisa CI. Sci.. 

[25] D. lonescu-Kruse: Elliptic and hyperelliptic functions describing the particle motion beneath 
small-amplitude water waves with constant vorticity, preprint. 

[26] D. Kinderlehrer, L. Nirenberg, and J. Spruck: Regularity in elliptic free boundary value prob- 
lems I, J. Anal. Math. 34 (1978), 86-119. 

]27] J. Lighthill: Waves in fluids, Cambridge University Press, Cambridge, 1978. 

]28] A. Lunardi: Analytic Semigroups and Optimal Regularity in Parabolic Problems, Birkhauser, 
Basel, 1995. 

[29] A.-V. Matioc: On particle trajectories in linear water waves. Nonlinear Anal. Real World Appl. 

11 (5) (2010), 4275-4284. 
[30] B.-V. Matioc: Analyticity of the streamlines for periodic traveling water waves with bounded 

vorticity. Int. Math. Res. Not. 17 (2011), 3858-3871. 
]31] J. B. McLeod: The Stokes and Krasovskii conjectures for the wave of greatest height. Stud. 

Appl. Math. 98 (1997), 311-334. 
[32] P. I. Plotnikov: Proof of the Stokes conjecture in the theory of surface waves. Stud. Appl. Math. 

108 (2) (2002), 217-244. 

[33] J. F. Toland: Stokes waves, Topol. Methods Nonlinear Anal. 7 [8] (1996 ]1997]), 1-48 ]412-414]. 
[34] E. Wahlen: Steady water waves with a critical layer, J. Differential Equations 246 (2009), 
2468-2483. 



I I 



